Quantitative prediction of effective toughness at random heterogeneous interfaces 
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The propagation of an adhesive crack through an anisotropic heterogeneous interface is considered. 
Tuning the local toughness distribution function and spatial correlation is numerically shown to 
induce a transition between weak to strong pinning conditions. While the macroscopic effective 
toughness is given by the mean local toughness in case of weak pinning, a systematic toughness 
enhancement is observed for strong pinning (the critical point of the depinning transition). A self- 
consistent approximation is shown to account very accurately for this evolution, without any free 
parameter. 
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Introduction - While the intimate link between sur- 
face energy of a material, and the resistance to crack 
propagation, or toughness, has been elucidated in the 
ideal cleavage case (Griffith's criterion), it was soon un- 
derstood that the same concept of crack toughness could 
be used in the more common situation where dissipative 
processes in the bulk of the solid (yet in the vicinity of 
the crack surface) contributes to (or even dominates over) 
the thermodynamic surface energy. A basic justification 
of this extension is that the additional dissipated energy 
is confined in a boundary layer within a small distance, 
S, from the crack surface, the process zone. For large 
scale structures, as compared to 5, an effective surface 
energy — integrating all dissipative phenomena — is an 
operational concept. 

In this up-scaling justification of the relevance of linear 
fracture mechanics, it should be emphasized that as the 
reference to a true surface energy is lost, so is the prop- 
erty that only the crack surface area matters. When a 
crack propagates in a heterogeneous solid, similarly, this 
heterogeneity may trigger different dissipative phenom- 
ena, resulting in a toughness which cannot be reduced 
to the bare surface energy. This effect have been ex- 
hibited in recent experimental studies [H-Q showing how 
periodic modulations of interface properties affect crack 
propagation. The modeling proposed by Tvergaard and 
Hutchinson Q showed that the interplay of local plas- 
ticity and crack propagation could also produce such an 
effect. 

Crack propagation and depinning transition - The ef- 
fect of heterogeneities on crack propagation has been 
early discussed in the context of reinforcement by crack 
trapping with tough particles 0-01 an d the detailed 
mechanism of overcoming a single obstacle has been clar- 
ified. This scenario leads to a local instability that 
can also be exploited to trigger dissipative (e.g. vis- 
cous) mechanisms that will result in an enhanced effec- 
tive toughness. Although the individual mechanism can 



be very precisely modeled and described, the difficulty 
of considering crack propagation in a random hetero- 
geneous medium is that "collective" pinning/depinning 
takes place and the wide distribution of avalanches (i.e. 
instabilities now involving potentially a large number of 
pinning sites) precludes a quantitative modeling of the 
effective toughness without using a detailed microscopic 
modeling. 

Indeed, the advance of a planar crack front through a 
disordered medium can be seen as a dynamic phase tran- 
sition The elastic energy release rate (or cquiv- 
alcntly the stress intensity factor) plays here the role of 
a control parameter. Below threshold, the crack front 
only advances over a finite distance and eventually stops 
(static phase). Above threshold, the crack front cannot 
be arrested any longer and it reaches a finite velocity 
(mobile phase). The effective (macroscopic) toughness 
corresponds to the critical threshold. Beyond the field of 
crack prop agation, the paradigm of the depinning tran- 
sition |16l | appears to describe as various physical phe- 
nomena as wetting of disordered interfaces, solid friction, 
plastic flow of crystalline or amorphous materials... The 
common characteristic feature of these different systems 
lies in the competition between the driving force, the 
pinning induced by a random potential at a microscopic 
scale and an elastic restoring force at large scale. 

Depending on the relative strength of the random po- 
tential and elasticity of the front, one usually distin- 
guishes two generic situations as illustrated in FigJTJ 
- in weak disorder conditions, the depinning front is only 
slightly perturbed and smoothly advances as a whole, 
with modest velocity fluctuations. Viscous dissipation 
can indeed be turned to arbitrary low values in quasi- 
static conditions, and hence only the bare average surface 
energy will be relevant for the macroscopic toughness. 
Disorder plays only a role on the geometry of the crack 
front. A theoretical description of the crack front geom- 
etry can be described through a perturbation expansion 
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about the homogeneous case (see e.g. [I?} for a recent 
review). Such an approach, in spite of the sophisticated 
theoretical treatment involved, cannot account for an en- 
hanced toughening. 

- in strong pinning conditions, the front advances inter- 
mittently, through a series of localized micro-instabilities, 
or "avalanches". During these unstable events, viscous 
dissipation is at play (even for quasi-static loading condi- 
tions) and finally contributes to the global energy balance 
as an enhanced effective toughness. A thorough discus- 
sion of the transition between weak pinning and strong 
pinning can for instance be found in [l8j . 

It is therefore crucial to characterize and quantify the 
onset of strong pinning and to propose a way to evaluate 
toughness in this regime quantitatively without having 
to recourse to an exhaustive numerical modeling. Such 
is the purpose of the present study, which is based on 
a self-consistent approach that was proposed in 19|, 2{| , 
but never tested quantitatively. 
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FIG. 1. Snapshots of crack fronts for different toughness land- 
scape illustrating the transition from weak (left) to strong 
(right) pinning regimes - same distribution, different corre- 
lation lengths: £ x = 16 (left), £ x = 1 (right). Two sets of 
two fronts separated by a small time lag are represented. In 
weak pinning conditions (left), the distance between succes- 
sive fronts is nearly constant, indicating that the local velocity 
only weakly fluctuates. In strong pinning conditions, parts of 
the front remain stuck while other ones have significantly ad- 
vanced, indicating a jerky dynamics. 

Numerical model - We consider here an interfacial 
crack front propagating in the (z, x) plane in x— direction. 
The location of the crack front at time t is denoted h(z,t). 
A random toughness landscape K c (z,x) is considered. 
We denote its average K c and its standard deviation 
S c . In the following, the problem is made dimensionlcss 
by scaling stress intensity factors and toughness relative 
to K c . The reduced toughness (denoted as lowercase) 
k c (z,x) — K c (z,x)/K c is now a random field of unit 
mean and standard deviation a = S c /K c . It is assumed 
to be short-range correlated and its correlation lengths 
are denoted £ x and £ 2 in the direction of propagation 
and orthogonal to it respectively. 

The average stress intensity factor along the crack 
front, Ko(t), i.e. the driving force, or its dimension- 
less counterpart fco(i) = Ko(t)/K c , depends on both the 



macroscopic loading and problem geometry. We consider 
the classical case of a displacement controlled mechanical 
test driven at constant speed. To account for the stiffness 
of the system (specimen and loading device) , the bound- 
ary condition is described by a slow and steady loading 
rate such that on average the crack front velocity is set 
to vq. For this aim, we introduce the stiffness e such that 



ko(t) =e[v t-h(t)] 



(1) 



where h(t) is the average position of the front at time t. 

The roughness of the crack front gives rise to a mod- 
ulation of the local stress intensity factor. Using a first 
order expansion in the roughness Q the dimensionless 
stress intensity factor at point (z,x = h(z,t)) can be 
written k[z, h(z, t)] = ko(t) [1 + k e i (z, h(z, t))] where 



k e i{z,h{z)) = — 

Z7T 



h(z) - h(z') 
[z - z'Y 



dz' 



(2) 



(J stands for the principal value of the integral) . 

Neglecting inertial effects, we consider in the follow- 
ing an over-damped dynamics and study the behavior of 
the crack front at the verge of propagation from above 
i.e. wo — > + . In that framework, the front obeys the 
following equation of evolution: 

d t h(z, t)=TZ [k (t) + k e i(z, h{z, t)) - k c (z, h{z, t))] (3) 

where 1Z denotes the positive part, lZ(x) = max(a;,0). 

Note that the time scale has been set so that the vis- 
cosity is scaled to unity. Two parameters thus remain 
that characterize the driving dynamics: the (reduced) 
stiffness e and the (reduced) velocity vq. In the following 
only the quasi-static limit vq — > + is considered. The 
random landscape of (reduced) toughness of unit mean 
is characterized by the standard deviation a and the cor- 
relation lengths £ x and £ z . Note that, in the quasi-static 
limit, this simulation setup is expected to give the same 
results than extremal dynamics simulations [21 1. 
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FIG. 2. Left: Effective toughness as a function of the cor- 
relation length £ x for different type of disorders and system 
widths L z . Right: Toughness distribution along the crack 
front that propagates in an exponentially correlated Gaus- 
sian landscape of width L z = 64 for different £ x , ( a — 0.125). 
Symbols and lines correspond to simulation data and theoret- 
ical predictions respectively. 

In the following we generated random landscapes of 
sizes L z x L x with a unit lateral correlation length £ z = 1 
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(the discretization length scale in the lateral direction) 
and a tunable correlation length £ x in the direction of 
propagation. Three different types of random fields were 
considered: Uniform- Gaussian (UG), Gaussian- Gaussian 
(GG) and Gaussian-Exponential (GE) where the first 
term/letter refers to the probability distribution func- 
tion of fc c and the second one qualifies its pair correlation 
function. The UG and GG disorders consist of grids of 
N z (= Lz/t^z) x N x (= L x /£ x ) random numbers from uni- 
form and Gaussian distribution respectively. The space 
between grids points in the x direction follows a uniform 
distribution of mean £ x and standard deviation ^ x /2^/3. 
The generation of the GE landscape consists of realiza- 
tions of an exponentially correlated Gaussian noise com- 
puted according to first order scheme due to Fox et al. Q 
with an integration step A = £z/50. In all cases, the 
continuous toughness landscape k c (z, x) is interpolated 
linearly between two grid points in the x direction. 

The standard deviation cr was varied in the range 
[0.125 — 1] and the correlation length £ x in the range 
[0.006 — 800]. We used a reduced spring constant e = 1 
and a reduced velocity vq = cr/ 20 (checked to be small 
enough not to influence significantly the results). Peri- 
odic boundary conditions along z are considered. Inte- 
gration of Eq. [3] is performed according to an explicit 
mid-point scheme. The time step St is chosen so that the 
maximum front increment is less than on tenth of the 
noise discretization size in the x direction. Starting from 
a flat configuration, we first integrate the crack front dy- 
namic over a de-correlation length equal to £ x (L z /t; z ) - 5 
in order to reach a statistical steady state. The effec- 
tive toughness k e s is then computed as the mean value 
of the driving force minus the driving velocity fco — vq 
along a propagation length equal to L x = 1024 £ x . The 
effective toughness k e g is finally averaged again over ten 
simulations (different statistical samples). 

Numerical results - The weak-to-strong pinning tran- 
sition induced by the shortening of the toughness cor- 
relation length has an important consequence: a spec- 
tacular increase of the effective toughness (the depin- 
ning threshold). Fig. [2](lcft) shows the dependence of 
the average value of the reduced stress intensity factor 
k e s measured along the propagation of the front for dif- 
ferent disorder types and the correlation length £ x of the 
anisotropic toughness landscape. For large values £ x we 
obtain (k c g) w 1 i.e. the mean of the toughness disor- 
der {k c }- However for low values of £ x , we observe that 
k e g departs from (k c ) = 1 and can reach significantly 
higher values. We also note that the value of the effec- 
tive toughness clearly depends on the type of distribu- 
tion and correlation of the microscopic disorder. Such 
a dependence is certainly not surprising since the effec- 
tive toughness is interpreted here as the threshold of a 
(dynamic) phase transition and as such is expected to 
depend on the microscopic details. We also note a slight 
but clear dependence on the system size: the larger the 



system, the larger the effective toughness. The enhance- 
ment of effective adherence is also reflected by changes in 
the toughness distribution (weighted by time) visited by 
the crack front as illustrated in Fig. [2] The full toughness 
distribution is clearly biased towards higher toughness 
values as £ x decreases. 





















X 






Lz 






30 40 50 



FIG. 3. Left: a band of random toughness in a homogeneous 
landscape traps the crack front; the lower the line stiffness, 
the larger the front deflection. Right: Graphical solutions of 
the implicit equation Afc e ; (x) — Ak c (x) which determines the 
equilibrium position of the crack front. The stable positions 
in the band are highlighted (in red and green) for two values 
of the line stiffness. The slope of the fine lines between the 
stable portions corresponds to the stiffness S of the crack 
front (a) For large line stiffness, the trajectory corresponds 
almost to the toughness in the band (weak pinning), (b) For 
low stiffness, the equilibrium position toughness distribution 
is skewed toward higher toughness (strong pinning). 



Interpretation: a self- consistent approach - This evo- 
lution can be rationalized in the framework of a self- 
consistent approximation proposed in Ref. [l9[ for illus- 
tration purpose, and never checked quantitatively. The 
top panel of Fig. [3] shows a toughness landscape con- 
sisting of a unique band of fluctuating toughness k c {x) 
in an otherwise homogeneous medium of toughness fco. 
Because of the elastic coupling, the crack fronts un- 
dergoes a deflection Aft = h — h(Q) proportional to 
the toughness contrast Afc = k c — kg. The associ- 
ated line stiffness S = Afc/ A ft can be analytically com- 
puted from Rice formula [l^, and is shown to scale as 
S oc l/£ z log(L z /£ z ). An effective medium approxima- 
tion, in the spirit of the self-consistent approximation, 
consists in evaluating fc e jj as equal to the value of fco such 
that the average deflection of the front taken over the en- 
semble of successive stable positions of the front is zero. 

From weak to strong pinning - The simplicity of the 
one dimensional picture of Fig. [3J allows one to define a 
criterion that determines the pinning regime. In absence 
of driving force, an equilibrium configuration of the front 
at position h(z = 0, t) = x in the band is obtained when 
the elastic restoring force, Ak e i(x) = SAh(x), balances 
the toughness contrast Afc c (cc). Depending on the respec- 
tive amplitude of the line stiffness and the local toughness 
gradient, one can obtain for this implicit equation either 
a unique solution or multiple solutions for the front po- 
sition x. A simple criterion for the transition from weak 
to strong pinning can thus be drawn from the onset of 
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multistability. Under these conditions, strong pinning is 
obtained when 



> 1. 



(4) 



As schematically shown in Fig. Etright) different sets 
of stable positions along the band are generated as a 
function of the line stiffness. For large S£ x /a the sta- 
ble trajectory follows closely the toughness of the band 
as expected from the weak pinning regime. Conversely 
for lower S^ x /a the crack front will only visit a subset of 
high toughness values that characterize the strong pin- 
ning regime. As a consequence, the distribution of tough- 
ness at stable positions is skewed toward higher values. 

Direct numerical integration of the effective medium 
model was performed in order to test quantitatively its 
reliability. S is first determined numerically to account 
for the discreteness of the simulated crack front. The re- 
duced toughness distributions and their means are then 
generated from one dimensional trajectories as exempli- 
fied in Fig. GUright). The comparison between the self- 
consistent approximation and the crack front simulations 
reported in Fig.[5]and in Fig.|3]shows a remarkable agree- 
ment. k e g is very accurately reproduced as a function of 
£ x , L z . a and the different disorder types. Note that 
the model not only accounts for the effective toughness 
variations but also for the visited toughness distributions 
without any free parameter. 
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FIG. 4. Effective toughness as a function of the correlation 
length £ x for different toughness standard deviations a and 
system widths L z . The larger a and/or the shorter £ x , the 
larger the effective toughness. Symbols and line correspond 
to simulation data and theoretical predictions respectively. 

Scaling analysis - One can show that only a combi- 
nation of scaled parameters really contributes to the ob- 
served results. If the toughness k, and distance x along 
the propagation direction, are transformed into 



k = (k-(k c ))/a 

X - x/^x 



(5) 



then the elasticity of the crack front line is S = S^ x /a. 
Note that k has now a zero average and unit variance, 
and that the criterion for weak to strong pinning is sim- 
ply S = 1. Fig. H] shows the reduced effective toughness 
k against the depinning control parameter 1/5 = a/ S^ x . 
All numerical estimates of effective toughness can be col- 
lapsed onto a single master curve that captures the tran- 
sition. For small values of the control parameter 1/iS <C 1 
the weak pinning regime holds, fc e ff = T whereas, for 
large values 1/S> 1, fc c ff significantly increases (strong 
pinning). Note that the system size dependence is here 
captured in the line stiffness parameter (the latter in- 
creases as the logarithm of the crack front length) . This 
has no effect in the weak-pinning regime, but introduces a 
rather counter-intuitive result in the framework of frac- 
ture mechanics: the larger the system, the smaller the 
line stiffness and hence the larger the effective toughness. 

It thus appears that the disorder (local toughness pdf 
and spatial correlation) has a dramatic effect on the pin- 
ning conditions. A crack front can encounter weak or 
strong pinning depending on its propagation direction. 
In the same spirit, if the toughness landscape is non- 
symmetrical, the effective toughness in a given direction 
will not be the same as in the opposite one. Our ap- 
proach thus offers a natural interpretation to the recent 
results presented in Rcf. [§], but also a simple theoreti- 
cal framework from which a detailed optimization can be 
carried out. 

Conclusion - The depinning of a long-ranged elas- 
tic line appears to be strongly dependent on the spatial 
correlation of the disordered landscape through which 
it propagates. The transition between weak and strong 
pinning is well captured by a simple criterion build on 
the toughness gradient in the direction of propagation 
and the line stiffness. A simple self-consistent approxi- 
mation describes very accurately the progressive depar- 
ture of the depinning threshold from its weak pinning 
value (the mean value of the disorder) to the higher val- 
ues measured in strong pinning conditions. In particular 
this approximation appears to accurately capture the de- 
pendence of the depinning threshold to the finer micro- 
scopic details (statistical distribution of toughness and 
spatial correlation function) . Let us emphasize that such 
a result may not have been expected from the fact that 
strong pinning involves micro-instabilities and collective 
phenomena (avalanches, ...) due to the underlying dy- 
namical phase transition. 

Bearing in mind that the enhancement of the effec- 
tive toughness does not originate from the initial land- 
scape toughness distribution only but also from its spa- 
tial correlation, our results open a promising way to de- 
sign anisotropic and tough interfaces. The present model 
gives also significant insights for safety calculation since 
it reproduces accurately the toughness distribution from 
which extreme statistics can be deduced. Finally, our ap- 
proach offers interesting theoretical perspectives for other 
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systems presenting the prototypical situation where elas- 
tic interface are trapped by random disorder. 

S.P. acknowledges the support of ANR project 
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